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In this work we study the effects of strong magnetic fields on hybrid stars by using a full general- 
relativity approach, solving the coupled Maxwell-Einstein equation in a self-consistent way. The 
magnetic field is assumed to be axi-symmetric and poloidal. We take into consideration the 
anisotropy of the energy-momentum tensor due to the magnetic field, magnetic field effects on 
equation of state, the interaction between matter and the magnetic field (magnetization), and the 
anomalous magnetic moment of the hadrons. The equation of state used is an extended hadronic 
and quark SU(3) non-linear realization of the sigma model that describes magnetized hybrid stars 
containing nucleons, hyperons and quarks. According to our results, the effects of the magnetization 
and the magnetic field on the EoS do not play an important role on global properties of these stars. 
On the other hand, the magnetic field causes the central density in these objects to be reduced, 
inducing major changes in the populated degrees of freedom and, potentially, converting a hybrid 
star into a hadronic star. 


I. INTRODUCTION 

Neutron stars undoubtedly belong to the most suitable 
environments for studying properties of strongly inter¬ 
acting matter under extreme conditions. For example, 
the density inside these objects can reach values much 
higher than the nuclear saturation density ~ 2.7 x 10^'’’ 
g/cm^. This makes neutron stars natural laboratories 
where one can examine and shed some light on the still 
open question concerning the equation of state (EoS) for 
ultra-dense matter and the role played by exotic degrees 
of freedom, such as hyperons and quarks. 

Another important feature present in compact objects 
and studied in this work is their strong magnetic field. 
From pulsar observations, the magnitude of the surface 
magnetic field in neutron stars has been found to be gen¬ 
erally of the order of 10^^ — 10^^ G. However, according 
to observations of star periods and period derivatives, 
classes of neutron stars known as Anomalous X-ray Pul¬ 
sars and Soft-Gamma-Ray-Repeaters can have surface 
magnetic field as large as 10^^ — 10^® G. These are usu¬ 
ally referred to as magnetars [IHi . One expects to find 
even stronger magnetic helds inside these stars as already 
calculated in Ref. According to virial theorem argu¬ 
ments, which give an upper estimate for the magnetic 
inside neutron stars, they can possess central magnetic 
fields as large as IOI 8-20 q 

The origin of strong magnetic fields in compact stars 
is still unclear. One common hypothesis involves the flux 
conservation of the progenitor magnetic field [HI- How- 
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ever, this idea is not suitable for magnetars since a canon¬ 
ical neutron star M ~ 1.4 Mq would require a radius less 
than its Schwarzschild radius in order to generate a sur¬ 
face magnetic field of the order of lO^^G [l^. Another 

S issibility suggested by Thompson and Duncan in Ref. 

describes a newly born neutron star combining con¬ 
vection and differential rotation to generate a dynamo 
process which is able to generate fields as large as 10^^ G. 
However, this standard explanation fails when trying to 
explain the supernova remnants associated with these ob¬ 
jects [ 1 ^ . 

Whatever the origin of strong magnetic fields might 
be, they affect locally the microphysics of the equation of 
state (EoS), as for example, due to the Landau quantiza¬ 
tion of the energy levels of charged particles and the effect 
of the anomalous magnetic moment (AMM) of charged 
and uncharged particles. Globally, magnetic fields affect 
the structure of neutron stars through the Lorentz force 
associated with the macroscopic currents that generate 
the field. They also affect the structure of the space- 
time, as magnetic fields are now a source for the gravita¬ 
tional field through the Maxwell energy-momentum ten¬ 
sor. As a consequence, magnetized stars are anisotropic 
and require a general-relativity treatment beyond the so¬ 
lution of the Tolman-Oppenheimer-Volkoff (TOV) equa¬ 
tions Bd. 

In this work, we model magnetized hybrid stars in 
a self-consistent way by solving Einstein-Maxwell equa¬ 
tions in the same way as done in Refs. d d- In 
Ref. d: the authors studied magnetized quark stars 
taking into account the magnetic field in the equation 
of state, the magnetization term for the matter, and also 
the magnetic field in the gravitational field equations. 
They have found that neither the magnetic field nor the 
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magnetization change significantly the global properties 
of these stars for a magnetic field strength of the order of 
^ 10^® G. Note that, the equation of state used by those 
authors, namely, a quark CFL model [l^ possesses a high 
baryon density range reaching two times the saturation 
density at the surface of star. 

In order to assess the role that a magnetic field de¬ 
pendent equation of state and the magnetization play in 
the global properties of stars, we use in this work a more 
complex equation of state than in Ref. [i3- Our EoS 
describes magnetized hybrid stars containing nucleons, 
hyperons and quarks and takes into account the anoma¬ 
lous magnetic moment for all hadrons. As a consequence, 
it produces a magnetization much higher than the one 
used in the Ref. [l^- Despite this, we show that the neu¬ 
tron star structure, like its mass-radius relationship, is 
not modified drastically with the inclusion of the mag¬ 
netic field in the EoS and the magnetization. On the 
other hand, the particle population is significantly mod¬ 
ified when the magnetic field is included. The main im¬ 
pact observed is the conversion of a non-magnetized hy¬ 
brid star with hadron and quark degrees of freedom to a 
highly magnetized hadronic star composed simply by nu¬ 
cleons. In reality, the temporal star evolution goes in the 
other direction, as the magnetic field of the star decays 
over time allowing a hadronic star to become a hybrid 
one. 


II. FORMALISM 


In this work, we study stationary highly magnetized 
neutron stars. In this context, stationary neutron stars 
with no mag netic-field-dependent EoS were studied in 
Refs, [l^, [HI- An extension of these works, where the 
authors included magnetic fields effects in the EoS, was 
presented in Ref. [i3- Here, we follow the same general 
relativity formalism and setup as in Ref. m- Details 
of the relevant equations, numerical procedure and tests 
can be found in the references cited above. We present 
here only the key equations that are solved numerically 
for the sake of completeness and better understanding for 
the reader. 

The choice of the coordinates in general relativity is 
crucial not only to write the gravitational equations in an 
advantageous form but also to make them easier to solve 
numerically. In the present case, due to the symmetry 
of the system, a polar-spherical type coordinate system 
is chosen, namely, the Maximal-Slicing-Quasi-Isotropic 
coordinates (MSQI) (for a review see Ref. The 

metric in the MSQI coordinate system is written as: 

ds^ = - sin^ 6l(# - N'^dtf 

- A'^idr'^ +r'^de'^), (1) 

with N^{r,9) being the shift vector, N{r,0) the lapse 
function and A and B functions of r and 6. The energy- 


momentum tensor of the system reads: 

={e+p)u^^u’' 

777 

+ -{b>^b’^ + gn) 

+ — + (6 • • 6 )) , ( 2 ) 

Mo V ^ / 

where m and B are the lengths of the magnetization 
and magnetic field 4-vectors. Here, the first and sec¬ 
ond terms correspond to the matter contribution. They 
were written separating the isotropic matter contribution 
(first term) from the magnetization contribution which 
is anisotropic (second term). The third term in Eq. ([2]) 
is the pure electromagnetic contribution to the energy- 
momentum tensor. The energy density is e, the isotropic 
contribution to the pressure is p and the fluid veloc¬ 
ity is The magnetic field in the fluid rest frame 
is b^ = diag (0,0,0, B), where B is defined to point 
into the z-direction. The metric follows the convention 
g^'^ = diag (—1,1,1,1). In the rest frame of the fluid 
(using = diag (I, 0,0,0)) we can write: 


'^mo — c, 

=p-mB, 

Tm2 = p-mB, 

T^3=P^ (3) 


where the components ^ and 2 9-^® usually referred 
to as perpendicular or transversal pressure (with respect 
to B) and as parallel or longitudinal pressure, re¬ 
spectively. The third term in Eq. ([2]) corresponds to the 
pure magnetic field contribution. It can be written as: 


TO =2^ 
2Mo’ 

1 

2po ’ 

52 
2mo 


rTi2 _ 

^B2 — 


rp3 _ 


(4) 


Following the notation of Ref. , the standard stress- 
energy tensor for the perfect fluid (PF) gives an energy 
density: 


EiPF) ^r2(e-hp)-; 


(5) 


while the momentum density flux is given by: 

=T\e + p)U. (6) 

The 3-tensor stress components can be written as: 

S(P/)r = =p, 


( 7 ) 
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and 

S^''P^=P+{E(P^^+p)U\ ( 8 ) 

where the Lorentz factor is given byF = (1 — [/^)“^ 
being U the fluid velocity defined as: 


u = ^(n-N*). (9) 

with the lapse function N‘^ and the angular velocity n, 
as measured by an observer at infinity (see Ref. for 
more details). The expressions for electric and magnetic 
field components as measured by the observer (Oq) with 
velocity can be written as 


E 


Oi 


Fapn^ 



or 


dA^ 

dr 


1 

N 


— + N* 

89 


dA, 

86 



being r^q the coordinate equatorial radius. The total 
gravitational mass of the star M is given by: 


M = J A^Br^x 

[N{E + 5) + 2N'^B{E + p)Ur sin 9] sin 9drd9d(j), (18) 


with E being the total energy density of the system E = 
and S'! total stress S'* = + 

where i = r, 9 and ij). 

Assuming that the matter inside the star has infinite 
conductivity, the electric held measured by the coming 
observer Oi (observer with velocity u^) must be zero, i.e. 
E'^ = u^Fajs = 0. This condition leads to the following 
relation of the 4-vector potential inside the star: 


dAt _ ^ dA^ 

dx^ 8x^ 


(19) 


which for the case of rigid rotation, fl = const, one has: 


Ea — 2 


= 0 


1 8A. 


1 8A^, 


’ Br"^ sin 9 89 ’ B sin 9 dr 


,0 


( 11 ) 


The corresponding energy density (electromagnetic 
magnetization) reads: 


E{em) = 1 r/j ^ 2^) E^E^ + B^B, 

2^0 -o/ 


( 12 ) 


At = -flA^. (20) 

As calculated in Ref. [l^, the equation of motion 
= 0) reads 


ox^ e-\-p 


^f) 


8A^ 

dx'- 


0 , ( 21 ) 


with V = v{r,9) := ln{N). The relativistic log-enthalpy 
H{r,9) in Eq. ([21]) is written as: 


and the momentum density flux can be written as: 


JEM) 

"'0 


1 

do 


A^{B^E^ 


E'-B^) + "^B^B.U 

ID 


(13) 


The 3-tensor stress components are given by: 


g{EM)r ^ 1 

2/ro 

2mB'^Bg 


B^Br) + 


(14) 


MEM)e 
^ 8 


= — {E'^Er - E^Eg + B^Br 

2pA 

2m B^Br 

“b f^’ 


B^Bg) + 


(15) 


H := In 


e+P \ 

mbUbC^ J ’ 


( 22 ) 


where mj, denotes the mean baryon mass 1.66 x 10“^^ kg 
and Ufi the baryon number density, respectively. As 
shown in Ref. [T^, one can relate the electric current to 
the electromagnetic potential through the current func¬ 
tion /: 


f - nf = {e+p)f [AA) , (23) 

and by integration of the equation of motion = 

0 ), one obtains: 

H (r, 9) + V (r, 9) — InT (r, 9) + M (r, 9) = const, (24) 

being M(r,9) the electromagnetic term induced by the 
Lorentz force 0 written as: 


q{EM)cI> 


(16) 


Besides the above definitions, other important 
quantities measured by Eulerian obeserver are the 
circunferencial radius Rcirc: 

Ecirc — E{r^q, ~pr^q^ ( 1 ^) 


M {r,9) = M {A^{r,9)) ■.= - [ f{x)dx. (25) 

JA,l,{r,e) 

The model built in this work used a current function 
/ = /o = const. As shown in Ref. [I^, other choices for 
/ are possible, however, they do not alter the conclusions 
qualitatively. For different values of fo, the electric cur¬ 
rent changes (see Eq. (1^ 1 and, therefore, the intensity 
of the magnetic field in the star changes. 
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III. MAGNETIZED EQUATION OF STATE 

When, in addition to a usual hadronic phase, neutron 
star models also include a quark phase, they are normally 
described by two different equations of state. These equa¬ 
tions of state are connected at the point in which the pres¬ 
sure of the quark phase becomes higher than the pressure 
of the hadronic phase. In our approach, we have instead a 
combined model of hadrons and quarks and one equation 
of state for both phases. In this case, we can study impor¬ 
tant features of the deconfinement phase transition, like 
the strength of the transition, the mixing of phases and 
also the accompanying chiral symmetry restoration. For 
this purpose, we extended the hadronic SU(3) non-linear 
realization of the sigma model [22l - [^ to include quark 
degrees of freedom in a spirit similar to the PNJL model 
in the sense that it uses the Polyakov loop <I> as 
the order parameter for deconfinement. The Lagrangian 
density of the model in mean field approximation reads: 

^ — ^mag + I^Kin + C,Int + i^Self + C,SB — U, (26) 

where besides the kinetic energy term for hadrons, 
quarks, and leptons (included to insure charge neutral¬ 
ity), the terms: 

jC-mag = “ Si V'i 

(27) 

Cint = - Sj + m<t>^ + 

(28) 

Fseif = -p mlp^ + 

+ 34 (cc4 + ^ + + 4^ + 2 ^) 

- A:o(cr2 -I- -I- (5^) - -f -l- 

- fc2 (4 + ^ + 3(72^2 + ^4^ _ _ J2)^ 

-fc4 ( 29 ) 

CsB = -mlUa - {yimlJk - C, 

(30) 

U = (oqT^ -I- -I- a2T‘^p^)^^ 

-ka3ro‘log(I-6$2-p8$3_3$4)^ (31) 

represent the magnetic and anomalous magnetic mo¬ 
ment (AMM) interactions with the fermions, the inter¬ 
actions between baryons or quarks and vector and scalar 
mesons, the self interactions of scalar and vector mesons, 
an explicit chiral symmetry breaking term (responsible 
for producing the masses of the pseudo-scalar mesons), 
and a potential U for the $ field. The later is impor¬ 
tant in order to reproduce a realistic structure for the 
QCD phase diagram over the whole range of chemical 


potentials and temperatures, including realistic thermo¬ 
dynamic behaviour at vanishing chemical potential as 
shown in Ref. [^ . 

In Eq. (I?7l) . qi is the electric charge of each particle 
in multiples of the electron charge e, is the elec¬ 
tromagnetic field, K represents the tensorial coupling 
strength of baryons with the electromagnetic field ten¬ 
sor, and = 1 ( 7 '^, 7 “^]/2 . The mesons included are 
the vector-isoscalars ui and cj) (vector meson with hidden 
strangeness), the vector-isovector p, the scalar-isoscalars 
a and C (scalar meson with hidden strangeness) and the 
scalar-isovector <5, with T 3 being twice the isospin pro¬ 
jection of each particle. The isovector mesons affect 
isospin-asymmetric matter and are, consequently, impor¬ 
tant for neutron star physics. The coupling constants of 
the model can be found in Ref. [l^, [Ig • The hadronic 
sector was fitted to reproduce the vacuum masses of the 
baryons and mesons, nuclear saturation properties, rea¬ 
sonable values for the hyperon potentials and the pion 
and kaon decay constants (/^ and fk). The quark sector 
was fitted to reproduce lattice QCD data at vanishing 
chemical potential and phase diagram information, such 
as the location of the critical end-point and a continuous 
first-order phase transition line that terminates on the 
zero temperature axis at around four times saturation 
density. 

The effective masses of the baryons and quarks are 
generated by the scalar mesons except for a small explicit 
mass term Mq and the term containing $: 

Mb = 9Ba<J + PbstsS + gscC + ^Ob + 9b<s>^'^, (32) 
M* = gqacr -k 9qST3S + 3gcC + + 3g$(l - $). 

(33) 

With the increase of temperature and/or density, the a 
field (non-strange chiral condensate) decreases in value, 
causing the effective masses of the particles to decrease 
towards chiral symmetry restoration. The field $ as¬ 
sumes non-zero values with the increase of tempera¬ 
ture/density and, due to its presence in the baryons ef¬ 
fective mass (Eq. (IS^ b suppresses their presence. On 
the other hand, the presence of the $ field in the effec¬ 
tive mass of the quarks, included with a negative sign 
(Eq. (l33l) '). ensures that they will not be present at low 
temperatures/densities. In this way, the interaction with 
the medium determines which are the degrees of freedom 
present in the system. 

The magnetic field in the z-direction forces the energy 
eigenstates in the x and y directions of the charged par¬ 
ticles to be quantized into Landau levels ip: 


EL = y + (^^ML + 2,pL\B - , (34) 

where ki is the fermi momentum and Si the spin of each 
fermion. The last term comes from the anomalous mag¬ 
netic moment (AMM) of the particle that splits the en¬ 
ergy levels with respect to the alignment/anti-alignment 
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of the spin with the magnetic field. The AMM also mod¬ 
ifies the energy levels of the uncharged particles 


El = (35) 

The AMM constants Ki have values Kp = 1.79, = 

— 1.91, KA = —0.61, = 1.67, Kj; = 1.61, = —0.38, 

K~ = —1.25, K~ = 0.06. The sign of Ki determines the 
preferred orientation of the spin with the magnetic field. 
For zero temperature, the sum over the Landau levels v 
runs up to a maximum value, beyond which the momen¬ 
tum of the particles in the z-direction would be imaginary 


^max 


E* 2 + s^k^B - Mf 


2\qlB 


(36) 


We choose to include in our calculations the AMM 
effect for the hadrons only, since the coupling strength 
of the particles Ki depends on the corresponding mag¬ 
netic moment, that up to now is not fully understood for 
the quarks. Furthermore, it is stated in Ref. that 
quarks in the constituent quark model have no anoma¬ 
lous magnetic moment, and in Ref. 2^, that the AMM 
of quarks from one-loop fermion self-energy is very small. 
For calculations including AMM effects for the quarks 
see Refs. [29l - [^ . The AMM for the electrons is also not 
taken into account as its effect is negligibly small. Prop¬ 
erties of the magnetized SU(3) non-linear realization of 
the sigma model were presented in Refs. [H, [s^ for an 
effective (ad hoc) variation of the magnetic field inside 
the star. 

For this work, the equation of state for charge neu¬ 
tral chemically-equilibrated matter was calculated at zero 
temperature and over a wide interval of densities and 
magnetic fields with very small steps. This 2-dimensional 
table was included in the Lorene C-I--I- class library for 
numerical relativity ©[ll in order for the code to find 
the correct magnetic field for each density in each direc¬ 
tion of the star and, then, find the corresponding values 
for the system’s thermodynamical properties. Our tabu¬ 
lated equation of state is available upon request. 


IV. RESULTS 


The equilibrium configurations are determined by the 
central enthalpy He and the choice of the magnetic dipole 
moment /r. We could have, instead, chosen a fixed cur¬ 
rent function /q and allow the magnetic dipole moment 
to vary. We chose the former in order to have a better 
control of the parameter space and to investigate exclu¬ 
sively the effect of the magnetic field on the star. In 
order to do so, we built equilibrium sequences for fixed 
magnetic dipole moments ^ defined as (see Ref. 0): 


2^cos9 . 

= ^(r) Ir-i-oo, 


(37) 


which is simply the radial component (the orthonormal 
one) of the magnetic field of a magnetic dipole seen by 
an observer at infinity. 


The magnetization is defined as M = —dfl/dB (for 
more details see Ref. [S^) and in Figured] we show: i) 
the case when the magnetic field is included only in the 
structure of the star (no EoS(B), no mag); ii) the effect of 
the magnetic field also into the equation of state on the 
neutron star structure without the magnetization term 
(EoS(B), no mag) and iii) the effect of the magnetic field 
also into the equation of state on the neutron star struc¬ 
ture plus the magnetization term (EoS(B), mag). We 
also show the non-magnetized cased denominated TOV. 



FIG. 1. Relation between the gravitational mass and the cen¬ 
tral enthalpy for non-magnetized and magnetized models. In 
the last case, we also include the effects of the magnetic field 
into the equation of state (EoS(B)) and the magnetization 
term (mag) . 


The results present in Ref. [13 show that the effects 
of the magnetic field in the equation of state and the 
magnetization on the neutron star structure are not con¬ 
siderable. Note, however, that those authors presented 
solutions for magnetized quark stars whose equation of 
state (CFL model) is very stiff and reaches, for instance, a 
number baryon density of about 2 times the nuclear satu¬ 
ration density at the surface of the star. In this work, we 
model not only the quark phase with up, down quarks, 
but also the hadronic phase containing the whole baryon 
octet, in a self-consistent way in the presence of the mag¬ 
netic field. 

The results in Figure [T] corroborate in part the results 
presented in Ref. [131 showing also that the magnetic field 
does not have a considerable effect on the maximum mass 
and radius of highly magnetized neutron stars through 
the effect on the equation of state for given magnetic 
fields. Only a very small reduction of the star mass is seen 
for the largest magnetic moment used in the most massive 
star calculated. We see, however, from our calculation a 
difference in the curves when the magnetization term is 
included. This is due to the fact that, in our case, the 
magnetization strength can reach a value of about 10 
times the value in Ref. [l3. Note that the effect of the 
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FIG. 2. Parallel (py) and perpendicular (px) pressure con¬ 
tributions (with respect to B) for different maximum gravi¬ 
tational star masses as shown in Figure [T] The perpendicu¬ 
lar pressure, px = P|| — mB, decreases as the magnetic dipole 
moment increases. 


magnetization is to decrease stellar masses. This fact is 
related to the negative sign in Eq. (ED. 

In Figure E] we show the EoS's associated with the star 
presented in Fig. [T] They are constructed following the 
magnetic field profile, which is generated assuming dif¬ 
ferent magnetic dipole moments. For the pressure in the 
direction of the magnetic field (parallel) the inclusion or 
not of the magnetization does not change visibly the mag¬ 
netic field profile and, therefore, the pressure in the star 
(the curves py (no mag) and py (mag) are overlapped). 
The parallel pressure does become larger (stiffer EoS) for 
larger magnetic moments, but such a fact does not repro¬ 
duce more massive stars. This stems from the fact that 
the magnetic field reproduced inside stars in this work is 
not large enough to turn the EOS much stiffer. In addi¬ 
tion, at the center of the star, where the magnetic field is 
larger, the baryon number density is reduced (for larger 
Bs) and the Haas-van Alphen oscillations make the EOS 
stiffer and softer at different densities [Sbl - fs^ . 

In contrast, the magnetization affects the equation of 
state contributing to the perpendicular pressure, which 
is reduced in accordance with the magnetic dipole mo¬ 
ment. Due to this, the maximum gravitation masses are 
also reduced. For example, the maximum gravitational 
mass for a dipole magnetic moment of /x = 3.5 x 10^^ Am^ 
is reduced from 2.49 Mq to 2.46 Mq when the magne¬ 
tization effects are included in the calculation. 

From this point on, all the results shown will in¬ 
clude, for consistency, the magnetic field effect in the 
EoS and the magnetization effect. The mass-radius di¬ 
agram for highly magnetized neutron stars determined 
by a constant magnetic dipole moment fi is presented 
in Figure El In this figure, we also show calculations 



FIG. 3. Mass-radius diagram for non-magnetized and mag¬ 
netized models. The calculation was done for different fixed 
magnetic moments /x. The higher the magnetic moment, the 
higher the magnetic field. Effects of the magnetic field into 
the equation of state and the magnetization are also included. 
The gray line shows an equilibrium sequence for a fixed baryon 
mass of 2.2Mq. The full purple circles represent a possible 
evolution from a highly magnetized neutron star to a non- 
magnetized and spherical star. 


for evolutionary sequences at fixed star baryonic mass of 
Mb = 2.2 Mq. These lines may represent the transition 
from a highly magnetized neutron star (a younger star) 
to a non-magnetized one (an older star). 

We have chosen a fixed baryon mass of 2.2 Mq be¬ 
cause its evolution line ends almost at the maximum 
mass for the non-magnetized and spherical configuration. 
Looking at different magnetic dipole moment lines, i.e 
/X = 1.0 X 10^^ Am^, 2.0 x 10^^ Am^, /x = 3.0 x 10^^ Am^ 
and /i = 3.5 x 10^^ Am^, one sees that increasing /x (and 
therefore the magnetic field) affects the structure of the 
neutron star in many ways. First, the maximum mass in¬ 
creases, but not so much as in the case of no-fixed baryon 
number, what had been already raised in Ref. using 
spherical approach. This is an effect of the Lorentz force 
acting outward and against gravity. For this reason, the 
star can support more mass. Second, the circular equa¬ 
torial radius of the sequence increases and the star be¬ 
comes much more deformed with respect to the symmetry 
axis. This deformation is also an effect of the assump¬ 
tion of a poloidal magnetic field, which makes the star 
more oblate. Calculations including toroidal magnetic 
field components have shown that magnetized stars be¬ 
come more prolate with respect to the non-magnetized 
case 

Figure E] shows the magnetic field profile and the en¬ 
thalpy iso-contours for a star with the maximum mass 
for a magnetization /x = 3.5 x 10^^ Am^, as can be seen 
in Figure El This value roughly corresponds to the solu¬ 
tion with maximum field configuration achieved with the 
code. 
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FIG. 4. Magnetic field surfaces (on the top panel), i.e 
iso-contours measured by the Eulerian observer Oo for the 
chiral EoS. This star is near the maximum equilibrium con¬ 
figuration achieved by the code and the maximum mass for 
the value /r = 3.5 x 10®^ Am^, as shown in Figure [3] In the 
bottom panel, the corresponding enthalpy profile is shown, 
which corresponds to a central enthalpy of He = 0.26 
(n = 0.463 fm~®). The gravitational mass obtained for the 
star is 2.46 Mq for a central magnetic field of 1.62 xlO^® G. 
The ratio between the magnetic pressure and the matter pres¬ 
sure in the center for this star is 0.793. 




0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

Bc(10’®G) 


0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

Bc(10’®G) 


1 

0.9 


■ / ■ 

J 0-8 

■ 


^ 0.7 

\ ■ 


0.6 

■ 

. 

0.5 

/ ■ 

..--r' . 

0.4 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
B,,(10'®G) 


0.1 0.2 0.3 0.4 0. 

Bc(10' 


5 0.6 0.7 ( 

■"s) 


FIG. 5. Microscopic and macroscopic star quantities, i.e 
central baryon number, gravitational mass, dipole magnetic 
moment and ratio between polar and equatorial coordinate 
radii, as a function of the central magnetic field. These 
curves represent an equilibrium sequence at fixed baryon mass 
Mb = 2.2 Mq for different magnetic field intensities. See the 
line for Mb = 2.2 Mq in Fig. [3]). 


Now we follow the fixed baryon mass Mb = 2.2 Mq in 
Figure [3] and present microscopic and macroscopic prop¬ 
erties as a function of the central magnetic field in the 
stars. As shown in Fig.[5l the central baryon number den¬ 
sity decreases with the magnetic field and it has the max¬ 
imum value at the center only in the static case. In the 
other cases, its maximum is found somewhere inside the 
star. This is due to the Lorentz force, which is related to 
the macroscopic currents that create the magnetic field, 
acting on the matter which has been pushed off-center. 
This is analogous to the number density reduction in the 
rotating star case. As we will see, this has a huge impact 
on the particle population of these objects. Second, the 
gravitational mass also increases, as already mentioned 
in this section. Third, as we fixed the baryon mass Mb, 
each star in the sequence reproduces different values of 
the central enthalpy He and a different magnetic dipole 
moment The latter is free to vary and increases with 
the magnetic field. Fourth, the ratio between the polar 
and the equatorial radii increases as the magnetic field 
increases and, therefore, the star becomes more deformed 
(oblate). 

For all quantities in Figure [Sj the curves have a quali¬ 
tatively change in behaviour for a magnetic field strength 
of 0.9—1.0 X 10^® G. At this point, the magnetic force has 
pushed the matter off-center and a topological change to 
a toroidal configuration can take place [10|. However, 
our current numerical tools do not enable us to handle 
toroidal configuration, which gives a limit for the mag¬ 
netic field strength that we can obtain within this ap¬ 
proach. As a consequence, the baryon number density, 
for example, will never reach zero at the center of the star. 
Still, the value of the magnetic field shown in Figure [5] 
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FIG. 6. Stellar particle population as a function of the baryon 
chemical potential. All figures represent an equilibrium se¬ 
quence at fixed baryon mass Mb = 2.2 Mq. As one increases 
the magnetic field, the particle population changes inside the 
star. These stars are represented in Fig. [3] by the full purple 
circles. For the non-magnetized case (B=0), the vertical red 
curve represents the chemical potential reached at the center 
of the star, namely, 1320 MeV. 


represents the limit in terms of magnetic field strength 
for a star at fixed baryon mass of Mb = 2.2Mq. Other 
configurations, as depicted in Figure 01 can reach values 
higher than 1.0 x 10^® G. 

Note that, in Figure [SJ the ratio between the polar 
and the equatorial radii can reach 50% for a magnetic 
field strength of ^ 1 x 10^® G at the center. Therefore, 
one sees that the deviations from spherical symmetry are 
quite significant and need to be taken into consideration 
while modelling these highly magnetized objects and a 
simple TOV solution can not be applied. 

The changes in the global properties of stars due to the 
inclusion of the magnetic field into the gravitational equa¬ 
tions are remarkable, and in order to study how the mi¬ 
crophysics is modified with the magnetic field, we present 
in Figure El the particle population Yi = pijpb as a func¬ 
tion of the baryon chemical potential pB for different val¬ 
ues of the magnetic dipole moment p for a fixed baryon 
mass Mb = 2.2 Mq . The kinks in the population plot 
that can be observed for values equal to or greater than 
p = 2.0 X 10®^ Am^ are due to the Landau quantization. 

For the spherical non-magnetized case, the TOV solu¬ 
tion was obtained for a hybrid star (with mixed phase) 
composed by the baryon octet, electrons, muons and u, 
d quarks (no s quark). In Figure El the red vertical line 
represents the baryon chemical potential reached at the 
center of the maximum mass star in the non-magnetized 
case. With the inclusion of the magnetic field through 
the dipole magnetic moment oi p = 1.0 x 10®^Am^, the 



FIG. 7. Equations of state for a magnetized and non- 
magnetized star with fixed baryon mass of Mb = 2.2 Mq. In 
the non-magnetized case, the equation of state describes a hy¬ 
brid star with mixed phase (without s quark). On the other 
hand, the magnetized star is described by an equation of state 
with neutron, protons, electrons and muons. 

central baryon chemical potential is reduced due to the 
Lorentz force. The new central value for /is is below 
threshold for the creation of quarks, which are, there¬ 
fore, suppressed. An even larger effect can be seen in 
the star for higher values of the magnetic dipole moment 
p = 2.0 X 10®^ Am^, p = 3.5 x 10®^ Am^, when even the 
hyperons are suppressed. As a result, the properties of 
these objects such as neutrino emission and consequently 
the star cooling, are strongly affected by the magnetic 
field strength in their interior as already pointed out in 
Ref. [13 for a spherical solution. 

In this way, younger stars that possess strong mag¬ 
netic fields might go through a phase transition later 
along their evolution, when their central densities in¬ 
crease enough for the hyperons and quarks to appear. 
Such phenomena might have observable signatures such 
as a distinct change in the cooling behaviour as well as 
a very different rotational slowing down, reflected in the 
stellar braking index. 

Figure [7] depicts the equations of state corresponding 
to the stars with B = 0 and p = 3.5 x 10®^ Am^ shown 
in Figure El For both cases, the equations of state are 
shown up to the maximum value of the energy density 
reached in center of the star. From Figure [71 the sup¬ 
pression of the mixed phase (which in the B = 0 case is 
in the range e = [3,4.5]fm“^) is evident. In the case 
p = 3.5 X 10®^ Am^, as already discussed, the Lorentz 
force reduces the energy density in center of the star to 
a value of ^ 0.8fm“^, i.e. 78% less that in the non- 
magnetized case. In this same figure, we show also con¬ 
tributions of the magnetic field and the magnetization 
to the equation of state. For the values of magnetic 
fields strength and conditions considered in this work, 
the former does not affect the equation of state. This 
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can be seen from the components of the parallel pres¬ 
sures, p|| (no mag) and py (mag), which are indistinguish¬ 
able. On the other hand, the magnetization has a small 
effect on the equation of state, which can be seen from 
the difference between the parallel and the perpendicular 
components of the pressures. 

V. CONCLUSION 

We studied highly magnetized neutron stars in a gen¬ 
eral relativistic framework. We solved Einstein-Maxwell 
equations self-consistently, including a poloidal magnetic 
field generated by macroscopic currents and taking into 
account all anisotropies associated with such a field. This 
includes the effects of the magnetic field on the equation 
of state, including the magnetization. In the EoS, we 
took into consideration the effects of the Landau quanti¬ 
zation of the charged particles in the presence of a mag¬ 
netic field and the anomalous magnetic moment of all 
hadrons (even the uncharged ones). 

We found that the leading contribution to the macro¬ 
scopic properties of stars, like mass and radius, comes 
from the pure field contribution of the energy-momentum 
tensor. Still, we found a small contribution when taking 
into account the magnetic field corrections in the equa¬ 
tion of state through the magnetization. In addition, the 
breaking of the spherical symmetry and the Lorentz force 
coming from the macroscopic currents, which generate 


the magnetic field, change the central density of stars. 
More precisely, high magnetic fields prevent the appear¬ 
ance of quark and mixed phases in highly magnetized 
stars. 

Assuming that the magnetic field decays over time, 
stars would not only become less massive and smaller 
over time, but also go through phase transitions to more 
exotic phases, thereby potentially significantly modifying 
their cooling behaviour, which will be studied in detail 
in future work. 

In the future, we are going to include efects due to 
star rotation in our calculations, which, for example, can 
effect the stellar braking index. Nevertheless, we expect 
that our findings are not going to change qualitatively, 
since rotation is going to further decrease the central den¬ 
sity of the star 
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